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DESCRIPTION AND TEXTS OF AUXILIARY PROGRAMS FOR PROCESSING 
VIDEO INFORMATION! PART II. SUODH PROGRAM OF AUTOMATED 
SEPARATION OF QUASIHOMOGENEOUS FORMATIONS 

By 

V. I, Borisenko and L. S. Chasalin 


Introduction / 3 * 

When ohservatlons are made £rom an aircraft or satellite the earth's sur- 
face seems to consist of individual sections, separated by various types of 
boundaries, roads, rivers, landmarks, tree -plantings, etc. 

In examining the color photographs, or those synthesized in conventional 
colors, of fragments of the earth's surface one can find sections that have the 
same color. Black-and-white photographs or quasiphotographs £l3 for each resolu- 
tion component yield an idea on the Intensity of the electromagnetic radiation 
in individual spectral zones reflected from the earth's surface. The system 
consisting of the recorder.' (human eye) and the analyzer (human brain), is 
capable of isolating individual section according to the color index, shape of 
the boundaries, texture, etc. This system is capable of grasping the entire 
picture as a whole or separating the image depending on the assigned or intui- 
tive criteria. The system ”eye-braln" can ignore individual sharp changes, and 
at the same time, note inconspicuous contours; it is capable of extrapolating 
open boundaries; it is easily adjusted to the different criteria of isolation 
or identification; it reacts differently to the changes in brightness or color 
index in individual sections of the imago (i.o,, it can exclude, or on the 
contrary, record the slightest changes). Knowledge of its operating principles 

^'^N umbers in margin indicate pagination in original foreign text. 
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would ptovid«» Invaluable assistance In constructing automated systems of video 
Infocmation processing, but since these principles are still not known, then 
different intuitive algorithms are of definite interest. This weak describes 
the program for separation of quasihomogeneous formations, l.e,, the formations /4 
whose spectral characteristics are constant on the average. 

Description of Algorithm and Technique for Separating Quasihomogeneous . 

Formations 

The Initial material for the operation of the given program is the video 
information in a standard color-superposition format [ 23 t where each resolution 
component on the earth's surface is described by a color vector (a*lP^ (i.j” 
number of lines and columns of imago component, k— number of spectral ranges in 
which the radiation rafloctad from the given fragment is recorded. 

In the definition of the quasihomoganoous formation [,33 the existence is 
implied of a coherent sot of imago elements whoso spectral characteristics 
differ little from the corresponding components of a certain vector (b°). In 
other words, the coherent sot of imago components is considered to be a quasi- 
homogenoous formation if all the points that depict each of its components in 
the measurement space are grouped faltly compactly, i,e,, occupy a small volume 
around the point (b®). 

Different methods are possible for determining the hypercenters (b®) of each 
quasihomogenoous formation: their a priori assignment, with tne help of compu- 

tation of the polymodal histogram peak centers, etc. We will dwell on the 
technique of a formalized determination of (b®) based on the hypothesis of the 
existence of formations that are quasihomogenoous in the spectral sign. 

It is evident from a definition of the quasihomofeenous formation [33 that 
in a certain line or column at least several neighboring components must exist 
with close spectral characteristics. After uniting such components into a group /5 
one can compute the mathematical expectation of their characteristics, that should 
be taken as (b ), Since the directions for an examination of the information 
according to lines or columns are not equivalent (since the video information is 
organized by lines on the magnetic tape), the following algorithm is suggested 
to deteamlne the initial group of elements belonging to a certain quasihomogenoous 
formation. After the next block (lino) of video information has been read, the 
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voctor amount of mathomatlcal axpactation is computed for the intensity of the first 
L components of the line, and the q.uasidispersion of the spectral charac- 
teristics of these components is determined i 



Hero k— amount of spectral canals, 1 — number of spectral canal, L-- quantity 
of components taken as the minimum possible in the lino for the existence of a 
quasihomoganaous forsiiation, j — beginning number of component in the lino of 

video information, k —weight coefficients, 

* 


After this the fulfillment is deter mined of the inequality that defines the 
degree of compactness in space of the multlzonal irtOasuramants of the selected 
sat of neighboring components in the line, i.a,, the correctness of the inequality 
is checked 

T> . ( 3 ) 


Here D is the maasuca of heteroganalty adopted for the given image and 
L 

corresponding to the average quasiseparation of the points from the hypercenter 
(b® ) . With fulfillment of inequality(3)the components from j to Jq+L are con- 
sidered to belong to a certain formation. If inequality (3) is not fulfilled 
the beginning number is Increased by a unit and the process is repeated until /6 
this inequality is fulfilled. The number of beginning components L (that par- 
ticipate in the computation b^ and D) will be called the number of dispersion 
components, or simply the dispersion length. 

After the inequality (3) Is fulfilled, the process begins of increasing the 
given homogeneous formation with the hypercenter (b )=(b ), Each subsequent 
component in the line is studied for its affiliation to the given formation 

with number sj the fulfillment of the following inequality is verified for it 

, ( 4 ) 
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where etc, , D — measure of heterogeneity, generally speaking, 

6 

not coinciding with D^, From our viewpoint It Is expedient for to be larger 
than Dt. The encountered component that does not satisfy (4) Is viewed as a 
random malfunction, while If there are two neighboring components that do not 
satisfy (4), the process of forming a quaslhomogeneous object with number s In 
the given line Is considered to be finished, and certain additional charac- 
teristics of the formation are defined. The table whose structure will be 
described below has Information on the isolated object, and the selection of 
a new quaslhomogeneous formation begins. For this j Is assigned a number that 
is a unit greater than the number of the last added component, and the process 
of searching for a new^ q.uaslhomogeneous formation is repeated again. 

The table that is formed after an examination of the next line, has the 
following appearances 



Its components have the following neaningj S — OjUantity/of components of 

s-th quasihomogeneous formation, N — its number in the line, CHAR and KD — 

® s s 

vector of the averaged spectral characteristics and vector of quasidispersion, /7 
j® , j® — beginning and final number of components of the given object, and 

i® — number of line in which the search for the homogeneous objects was made. 

The meaning of the values and j® will be explained below, 

s s 

After a review of the next line in the computer memory, analogous tables 
are simultaneously found for the current and previous lines (NEWTAB and OliCTAB). 

Now the quasihomogeneous formations for the column are united. For each com- 
ponent of the table OLDTAB a corresponding (for the column) component NEWTAB 
is sought, for which the fulfillment of the system of inequalities is verified: 




A/ew 


^ (^h) 040 


(5) 


i.e,, is verified if in (l+l)-th line the object that can become a continuation 
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of the ^-th ohjact of the i-th line. If there Is no such object, then the object 
with number Is considered to be complete, and information about it is trans- 
ferred to the table RESTiffl for subsequent classification, while if such an object 
does exist, than the identity of the objects 1 and s is verifledi for this the 
quasidistanca between their characteristics in the space of multlspectral measure- 
ments is computed i 

ft » ^ £ I I . (6) 

c*i 

If this distancQ is smaller th'rn the assigned measure of heterogeneity D , the 

c 

objects are united into one. This entails a certain reformation of the s-th 
component of NEWTABj 

1. The quantity of elements of the S formation is replaced by the sum S +S„| 

S S w 

2. The characteristics of the object CHARg and KDg recalculatad into average 
characteristics of the summary formationi for this their vector amounts are added 
with the weight coefficients that are proportional to the areas of the coherent 
objects f 

3. i® is replaced by , j® is assigned the value j^ from OLDTAB. (During the /8 
formation of NEWTAB the number of the extreme column of the Isolated object is 
entered into j®) . The presence of the numbers of linos i® and 1® , between which 

the given object is located, and the numbers of the extreme column of the upper 
line of the object j® facilitates the identification of the table with the 
initial material. Such coordinate correlation is necessary at the stage of 
further interpretation. 

4. The so-called "birth number" of the given object is entered into N^; it appears 

s 

during the formation of a new quasihomogenoous formation that is not a continu- 
ation of the object from the previous line. The presence of this number in NEWTAB 
indicates that this object is already linked to the objects of the previous lines, 
and after a review of the entire table OLDTAB, the next numbers are assigned to 
the formations that are not a continuation. 

The RESTAB table that is constructed analogously to NEWTAB or OLDTAB is formed 
from objects of OLDTAB that do not continue in the next line. Only j® and j® are 
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excluded from it, emd N is formed with a rewriting of the next component OLDTAB 

s 

into RESTAB. When the assigned number of components of RESTAB is filled, the 
information is released onto magnetic tape. Consequently, all the data on each 
object (axeai spectral characteristics! location on the imagei numbeis obtained by 
the given object during its emergence and release onto the magnetic tape) are 

preserved for the subsequent classification, 

• 

At the same time as the review of the next line, another output tape is 
formed that is necessary for the display of the isolated quasihomogeneous forma- 
tions. The length of the line of this tape is determined by the length of the line 
of initial mater iali two bytes are assigned to each color vector of the image i in 
these bytes the "birth” numbers of the given formation are entered. Thus, each 
line of the resulting magnetic tape is filled with numbers N^, from which one can /9 
always, after using RESTAB, determine the characteristics of each specific quasi- 
homogeneous formation. The absence of numbers in .certain components of the output 
tape indicates that these components do not belong to any of the isolated objects. 

After the end of the review of the fragment whose coordinates (number of 

beginning and final lines and columns) are assigned from the controlling punched 

cards, the total area is computed for the separated quasihomogeneous objects P, 

the percentage relationship is determined for this total area P to the area of 

the entire fragment S_ , and besides these two amounts, the number of separated 

objects # is also printed. This makes it possible (at the stage of determining 

the separation parameters) to evaluate the performance capacity of the given 

algorithm, and to make the appropriate corrections! increase or decrease and 

D , change L, etc. With the selected parameters the amount P shows the expediency 
s 

of automated Aeparationi the small value of P indicates either the incorrectness 
of the set task (search for quasihomogeneous formations where they do not exist), 
or that for the given fragment of the image, the parameters for separation ware 
not successfully selected (if it is known a priori that the given material contains 
fairly many homogeneous formations) . 

The volume of the operational memory that is occupied by the program to- 
gether with the regions of input-output Is 3*5 bytes, so that the algorithm can 
be used on a- computer with fairly small operational memory. The process of 
separation occurs on the ES-1020 computer at the rata of reading the magnetic tape 
with initial material. 
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Tha simplicity of ths algorithm, xaplacsmsnt of axpansiva (in tha sansa of 
fast-xasponsa) opaxations of multiplication and avolutlon (in tha formulas fox 
computing distancas) with simpla opaxations of tha moduli sum typo— all of this /lO 
dstorminos tha high spoad of tha program xoalizatlon of tha suggostod algorithm. 

As shown hy tha xosults of prollminary calculations, such a xoplacomant is q.ulto 
justlf;.ad, sines it providos a largo saving of axpondod machina tlma without 
strongly affecting tha accuracy of the separation. 

The SUODH program was written in tha ASSEKBIER language fox the ES computer t 
its block-diagram is given in the appendix, and its text in a separate appendix. 

It was computed to process video information with any number of spectral 
(no more than 8) and arbitrary line length, Tha size of tha RESTAB block is 
fixed (32 bytes), the number of blocks (objects) in one physical output tape 
recording is 32. 

2. Instructions for Working with tha SUODH Program 

The parameters necessary for operation of tha peogtam arai 
NST — beginning number of the line of the examined fragment. 

KST — final number of the line, 

NC3T,KCT — numbers of beginning and and of columns (color vectors), 

LCOL — number of spectral channels ’in studied material, 

Dj^—maasure of heterogeneity (under conditions of measurement units) adopted 
to form the homogeneous object, 

Dg — permissible deviation for the line necessary to determine tha affiliation 
of tha next component of the line to the formed homogeneous formation, 

Dq— measure of hatarogenaity that makes it possible to link tha objects of the 
neighboring lines into one object, 

LRAZ — number of dispersion components, 

Tha amounts of tha described parameters in tha form of four -value whole /ll 

numbers must be packed on the controlling card in the following order: 

L COL, MS T, MC T, LRf!2, ,7)i.KST, Kc r . 

Information input-output is written on a logical level [43, therefore, 
among tha controlling cards of tha assignment there must be perforated cards 

present that identify tha input -and output magnetic tapes, i.a,, operators 
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(f TTBL with indication of tho naraas of the input and output files. The program 
guarantees recording of the tail markers on "both output tapes after review of 
the last line of initial material indicated hy tho parameter KST. ' 

The numbers of the lines and the columns in the output tape correspond to 
tho numbers of the lines and tho columns of the Initial video information. 

Suggested Variants of Program Modification, 

1, The weight coefficients K® that serve to intensify the influence of 
individual channels participate in the algorithm for separation. In the des- 
cribed program variant their amount is tho same. However, a whols series of 
reasons (a priori information on the greater importance of one or several 
channels as compared to tho othorsi information on errors in tho operation of one 
of tho channels, etc) can result in tho need to change individual (or all) tho 
weight coefficients. Therefore, it is further suggested that weight coefficients 
bo introduced from a separate punched c&rd, 

2, In the operational computer memory tho tables NEMTAB and OLDTAB are assigned 

800 bytes each for placement of information on 25 objects that are located in 
one lino. However, in tho localities characterized .by groat diversity of natural 
formations, there can bo more homogeneous objects in one lino than Z5> I"t is /12 

suggested that these files be transferred to tho and of the program, and that the 

size of tho memory sat aside for them be expanded, 

3, In the examined algorithm for separation of tho q.uasihomoganoous formations, 
only the dispersion of the vector characteristics for the object was computed and 
stored in the memory. It is proposed that further perfection in the algorithm 
take into consideration tho moments of higher orders as w^ll. Their determination 
and storage req.uiras additional writing of tho computation blocki and increase 

in tho area of tho memory that is occupied by one component of tho table. 

Tho authors are deeply grateful to Z. A. Stotsonko for assistance in sotting 
up this work. 
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Appendix . Block Diagram of SUODH 



Beginning 

Identification of tapes, 
introduction of paiameters, 
establishment of 1-NST 


reading of next line of Infcrmation, 
sending of NEWTAB into OLDTAB, 
establishment of j"NCT 


dispersion block 
computation of D 


vD < Dt 


y . 

recSviry of next component of line 
computation of D 


storage of characteristics of 
I isolated ob^ct in MEWTAB 


unification of objects of neighboring 
lines, formation of lino of output 
magnetic tape, output of finished 
objects into RESTAB, release of RESTA; 
block to magnetic tape. 

- 


release of last RESTAB block to magnetic 
tape, output of results on alphabet-digital 
printer, entry of tail markers of output 
magnetic tape, i 
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